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Abstract 

We have generalized the curvature coupling models of Higgs inflation to study inflation with 
a scalar field coupling of the form ^f +2 ^_ 4 . Such higher order scalar and graviton terms can be 

expected from the renormalization close to Planck scale. We compute the amplitude and spectral 
index of curvature perturbations generated during inflation and fix the parameters of the model 
by comparing these with the Planck+WP data. We find that if the scalar self coupling A is in the 
range (10 -5 — 0.1), parameter a in the range (2.3 — 3.6) and b in the range (0.68 — 0.22) at the 
Planck scale, one can have a viable inflation model even for £ ~ 1. The tensor to scalar ratio r in 
this model is small and our model with scalar-curvature couplings is not ruled out by observational 
limits on r unlike the pure ^<3? 4 theory. By requiring the curvature coupling parameter to be of 
order unity, we have evaded the problem of unitarity violation in scalar-graviton scatterings which 
plague the Higgs inflation models. We conclude that the Higgs field may still be a good 

candidate for being the inflaton in the early universe if one considers higher dimensional curvature 
coupling. 
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1. INTRODUCTION 

The idea that the universe through a period of exponential expansion, called inflation 
[1-6] has proved useful for solving the horizon and flatness problems of standard cosmology 
and in addition providing an explanation for the scale invariant super-horizon perturbations 
which are responsible of generating the CMB anisotropies and formation of structures in 
the universe. A successful theory of inflation requires a flat potential where a scalar field 
acquires a slow-roll over a sufficiently long period to enable the universe to expand by at 
least 60 e-foldings during the period of inflation. There is a wide variety of particle physics 
models which can provide the slow roll scalar field 'inflaton' for inflation [7]. From the 
observations of CMB anisotropy spectrum by COBE and WMAP [8] it is not yet possible to 
pin down a specific particle physics model as the one responsible for inflation. In the light 
of recent discoveries by CMS [9] and ATLAS [10] it is of interest to consider the Standard 
Model Higgs boson as the candidate for inflaton. On the face of it the idea does not work as 
the inflaton quartic coupling should be of the order A ~ 10~ 12 to explain the amplitude of 
CMB perturbations measured by WMAP [8] while the 125 GeV Higgs has a quartic coupling 
A ~ 0.13 at the electroweak scale which can however go down to smaller values at the Planck 
scale due to renormalization [11-16]. However just from the standard model renormalization 
one cannot have the Higgs coupling A ~ 10~ 12 over the entire range of the rolling field (10-1 
)Mp during inflation and the standard slow roll inflation with a Higgs field does not give 
the observed amplitude and spectrum of density perturbations [17]. If the Higgs and top 
mass are fine tuned then there can be a small kink in the Higgs potential and the universe 
trapped in this false vacuum can undergo a period of inflation [18-20]. 

A way out of fine tuning the scalar self coupling to unnaturally small values was given by 
Fakir and Unruh [21] who pointed out that if one couples the scalar field to the Ricci scalar 
£$ 2 i? then the effective potential in the Einstein frame becomes a slow roll one with the 
effective scalar coupling being A/£ 2 . Density perturbations from inflation in the curvature 
coupled theories were calculated in [22, 23]. The equivalence of the density perturbation in 
Jordan and Einstein frame was shown by Komatsu and Futamase [24] who also calculated 
the tensor perturbations and showed that the tensor to scalar ratio is generically small in 
the Fakir- Unruh [21] model. 

Bezrukov and Shaposhnikov [25-27] revived the large curvature coupling model to mo- 
tivate the idea that the standard model Higgs field could serve as the inflaton in the early 
universe . The amplitude and spectral index of density perturbations observed by WMAP 
can be generated by the Higgs field with self coupling A ~ 0.1 and curvature coupling £ ~ 10 4 
[25-28]. This large value of £ needed however is seen as a problem as at the time of inflation 
the Higgs field is at the Planck scale and graviton-scalar scatterings due to the curvature 
coupling of the scalar would become non- unitary [29] . Ways of solving the unitarity violation 
problem in the Higgs inflation models have been explored in [30-42]. Higgs inflation with 
large scalar couplings have also been studied in [43, 44]. 

In this paper we assume that the dominant interaction between Higgs field and gravity 
is through operators of the form 

~ /Ma+26-2 
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The complete dynamics of the Higgs field involves the role of the Goldstone modes as has 
been studied in detail in [45-47]. The multifield dynamics of the Goldstone modes gives rise 
to sizable non-gaussianity. We will study the dynamics of the Higgs mode and impose a 
charge conservation and CP symmetry such that the Goldstone modes of the Higgs field do 
not acquire vevs. We will take the background Higgs field to be 

«=(£) p> 

where <£> is the Higgs mode with mass 126 GeV. Our inflation model falls in the class of 
inflation in /($, R) theories studied in Ref. [48]. Our motivation is that we use the Higgs 
quartic coupling A('H^'H) 2 where the standard model value of A(/x ~ M P ) can lie in the range 
A = (10~ 5 — 0.1) depending on the value of top quark mass [15, 16] or on new physics [49]. 
We take curvature coupling £ to be of order unity and check the possibility of generating 
the observed density perturbations from Higgs inflation by varying parameters a and b. 

We derive the curvature perturbation during inflation in two different ways. We de- 
rive the perturbations of modified Einsteins field equation in the Jordan frame in presence 
of the Higgs-curvature interaction terms and derive the amplitude and spectral index of 
curvature perturbation. We find that to generate the Planck+WP preferred amplitude 
= 2.1955+°J|j x 1(T 9 and spectral index n s = .9603 ± .0073 [50] for A = 10~ 3 we should 
have a ~ 3.02, b ~ 0.49 ( and for A = 0.1 we need a ~ 3.56, b ~ .22). In these fits we take 
£ = 1. Next as a consistency check we make a conformal transformation which removes the 
higher dimensional Higgs-curvature term but the effect of this term appears in the poten- 
tial in the Einstein frame. We calculate the curvature perturbation in the Einstein frame 
and find that our results are consistent with the Jordan frame results. Calculation of the 
curvature perturbation in both Einstein and Jordan frame for the £<3> 2 i? theory has been 
done previously in Ref. [51-54]. In the section (2) we derive the curvature perturbations 
and tensor perturbation in our theory in the Jordan frame and in the section (3) we make a 
conformal transformation to go to the Einstein frame and compute the curvature perturba- 
tions. Further in section (4) we discuss the renormaization of £$ 2 i? theory. Finally in the 
last section (5) we compare the results of the two different frames and discuss the viability 
of the Higgs inflation model in presence of higher order gravitational couplings. 



2. CALCULATION IN THE JORDAN FRAME 



In this section we introduce a scalar-gravity interaction term /(<&, R) in the action and 
calculate physical quantities related to the inflationary density perturbations such as the 
spectral index, curvature perturbation and tensor-to-scalar ratio. We start with the action 
for a scalar field interacting with gravity of the form 



S = J d 4 x^ 



where we take, 



f(Q,R) = R + 



2k 2 



/Lfa+26-2 
p 



A$ 4 



(3) 



(4) 
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where £ is a dimensionless coupling constant. Varying the action (3) w.r.t g^ v and $ we 
obtain the field equations, 

= FR^ - l -Jg^ - y^ u F + g^UF = k 2 (v^V^ - ^ V$V P $ - , (5) 

™ = - ^> (6) 

where F = Of /OR = 1 + fpg. 



2.1. Background quasi- De Sitter solution 

For the unperturbed background FRW metric (— 1, a 2 (t), a 2 (t), a 2 (t)) and scalar field $ = 
the above Eqs. (5) and (6) reduce to the form 

3FH 2 + l -(f-RF) + 3HF = k 2 Q</> 2 + (7) 



-2FH - F + HF = K 2 (f) 2 (8) 

L 

2k 1 

a nb — 1 



+ 3if0 + \/ -^ = . (9) 



Now we assume the second term of F i.e. is dominant for some values of a and b. 

We find this assumption to be valid while solving numerically for the values of a and b in 
our model which give rise to the experimentally observed density perturbations as shown in 
the Table (I). From Eq. (7), under this assumption and considering the slow roll parameters 
which are defined in Eq. (28) as small, the Hubble parameter in the Jordan frame turns out 
to be of the form 



4-q 
2b 



H= ^ r (±) M p . (10) 



From Eq. (9) under the slow roll assumption we get 



: ran . . 

* = "sit 1 " W^>1 < n > 



2.2. Scalar field and metric perturbations 

Now we perturb Eqs. (5) and (6) by perturbing the scalar field $ = 4>(i) + 5<j>(x, t) and 
the metric as 

ds 2 = -(1 + 2a)dt 2 - 2a(t)(d i f3)dtdx i + a 2 (t) (5^(1 + 2ijj) + 28^) dx l dx\ (12) 

where, a, if), (3 and 7 are scalar perturbations. We derive the Einstein equations for the 
f(R, 0) theory [55, 56] keeping the first order terms in the metric and scalar field perturba- 
tions. The component SGqq is given by 
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— — tb + HA = — 

a 2 {ty 2F 



3H i + 3H + jL^j 5F - 3H5F + l - - U) 5<f) 



{3HF - K 2 <p 2 )a + FA 



(13) 



and taking the difference 5G\ — 5G$ we get 



A + 2HA + [3H + 



A 

W) 



a 



2F 



35F + 3H5F — ( QH 2 + 



A 



8F + An' 



+ (-2kV + f t4> ) 8(f> - 3Fa - FA 

2 



(AK 2 <j) 2 



+ 3HF + 6F)a 



(14) 



where A = 3(Ha — tp) — A%/a 2 (t) and X — a (t)(/3 + aj). Here, in arriving the Eqs. (13) 
and (14), the leading order Eqs. (7) and (8) are also used. The other components SGio and 
SGij (i 7^ j) of the first order perturbed Einstein equation can be written as 



Ha — ip 



1 

2F 



and 



K 2 <j)5<p + 5F- H5F - Fa 



1 



(15) 



X + Hx-a-^j = ~[6F-F X ) 
respectively. The equation of motion of scalar perturbation is 



6<f>+3H5(/)+ 
where 



A { f2V 4> -f i<t) /K 2 



a 2 (t) 



<pa+(2<p + 3H^j a+j>A+^F t4 



(16) 



f ) ■ <-) 



5R=-2 



4 - 1 .4 H ( + 3H ] a + 2-^^P 



a?{t) 



a 2 (t) 



(18) 



Now we analyze the curvature perturbation 7Z — ip — H5<p/<p by choosing a gauge where 
5(f) = and <5i? = 0. This sets 7Z = ip and moreover we have SF = via 5F = [OF /dip) 8<p + 
(dF/dR) 5R. Under this gauge the Eq. (15) gives, 



n 



a 



H + F/{2F) 



(19) 
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and hence from Eq. (13) we get 



i / a ( ZHF - K 2 <j) 2 ) 1Z 

A = n _|_ _V / I C20) 

H + F/(2F) ya 2 (t) 2f(h + F/{2F) X ' ' 



Using Eq. (8) and Eq. (14), we obtain 

i (err F \ a 3F . f 3F + QHF + K 2 <j) 2 A \ 
A+ { 2H+ 2F) A+ 2F a+ { 2F - + W)) a = °- P1) 

Now we may write the differential equation for curvature perturbation by using the above 
Eqs. (19), (20) and (21) as 

•■ (a 3 (t)Q 8 )' k 2 



where, 



Q, = M (23) 



(H + F/(2F) 



In arriving Eq. (22), Eq. (8) is again used. Now one may re- write the Eq. (22) in terms of 
variables oj = a^fQ~ s and a k = ulZ as 



u" 



°k+ [k 2 - — )<T k = 0, (24) 



where prime denotes the derivative with respect to the conformal time defined as drj = 
dt/a(t) and 



OJ 



n „n 



a"(t) + <mo^ + i_Q:_ifq ± m 



oj a(t) a(t) Q s 2Q S 4 \Q S 

under quasi de-Sitter expansion a{rf) = gT? ^- £l ) anc ^ hence = ^ [2+3ei] and a'(t)/a(t) = 
a(t)H. Therefore we have 



u" 1 

where 



oj rf 



2 1 



(26) 



2 9 r 4 

V-T} — - 



3 

In arriving at the above expression we have defined 



1 + - (2ei + e 2 - e 3 + e 4 ) ] ■ (27) 



H (j) F E 

61 ~~ ~~tP ' 62 ~~ iJ0 ' 63 ~~ 2HF ' 64 ~~ 2HE ' (28) 
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E = F+ 3£l = QJ 1±I1 l 
2k 2 <P 2 <p 2 /(FH 2 ) 

Here are slow roll parameters and e'j terms have been neglected. Equation (24) then has 
solutions in the Hankel functions of order vr 



[d<>(*M) + 02H&>(k\r,\)] (30) 

Applying the Bunch-Davies boundary condition a(kr] — > — oo) = e tkv /\/2k we fix the inte- 
gration constants ci = 1 and c 2 = 0. Using the relation H v {k\q\) = -^r(i/) ^^jpj for the 

super-horizon modes hi] — > 0, we obtain the expression for the power spectrum for curvature 
perturbations is defined as 

^s^^y^r 1 (3i> 

The amplitude of the curvature power spectrum turns out to be 

1 f H 



IT (32) 



in 



and the spectral index is 

n n - 1 = 3 - 2v n ~ -4ei - 2e 2 + 2e 3 - 2e 4 ~ -6ei . (33) 

Using slow roll parameters, Eq. (23) can be simplified to the form Q s ~ QFe\M 2 with 

jr^2 << 6e| which will be justified later in the subsection (2.4). In our model of /($, R) 
coupling we find that to the leading order we have the relation e\ ~ — e 3 and e 2 ~ —e 4 
and these relations are used in the calculation of perturbation amplitude and spectral index. 
Plugging the values H and 4> from Eqs. (10) and (11) into Eq. (28), the slow roll parameters 
can be written as 



a + 2h-4 
6 



ei = b-\a - 4) (2 - ^-^(a + 26 - ^A^ 1 )/ 6 ^ (^J (34) 

a + 2b-4 

e 2 = r 1 (2-6) (1 - b)/b (a + 26-4)(a + 66-4)A (6 - 1)/b e 1/b (^ * ■ (35) 

For our model, we can write the expressions for the amplitude of power spectrum and the 
number of e-folding as 

A2 5[(2-5)/A]3-lMr^r^-^ m) 
n 288(a-4) 2 (a + 26-4)% 2 1 ' 

and 

(37) 



*' H b[{2-b)/\fr^ ( <j> 



Ni 

d> T 2(a + 26-4) 2 \M p 

4> f 

respectively. Here 0j and 0/ are the values of scalar field at the beginning and the end of 
inflation respectively. 
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2.3. Tensor Perturbations 



We define the perturbation of metric as follows 

9,u= 9,u + V and g^ = g^ + h^, 
where g^ u is background metric and 



h tJ = 



1 



a*(t) 



hij , h 



iO 



1 



hio, h 00 — —/loo- 



a?{t) 



(38) 



(39) 



Now to get the equation of tensor pertubation, we set h i0 = /i o = in the calculation. 
From the decomposition theorem, the non zero spatial components are traceless and 
divergenceless, i.e., 



ha = 0, dih i:j = 0. 
Using Eqs. (39) and (40), we obtain 

SR 00 = 0, 5R i0 = 0, 



(40) 



(41) 



SRij — 



1 2 , 1 - a • i a 

V h;j H — hij hij + 2 



2a 2 (t) ' 13 2 13 2a" 
So, perturbing Eq. (5), we obtain 

' lurDr, - I -lur - -aah f 



hi~, 5R = 0. 



(42) 



-Fa 2 A, 
2 3 



1 3 \ F 
-Fa 2 + -aaF V 2 D 

2 2 3 2 



2-F -2F 

. a 



a / a 2 



a 2 Aj, (43) 



where Aj = h^ja 2 . The right hand side of Eq. (43) vanishes by Eqs. (7) and (8). Thus we 
have 



(a 3 F) 



o. 



In the terms of polarization tensors ej- and e 2 -, the tensor Aj is written as 



A, = £>i4 + £> 2 e 



(44) 



(45) 



For gravity wave propagating in z direction, the components of polarization tensor are given 
by 



-e 1 =1 e 



2 =e 2 



1 1,2 1,2 n 



So the Eq. (44) can be written as 



A + 



(q 3 F) 
a 3 F 



A + — D A = 0. 



(46) 



(47) 
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where A = 1, 2. Now substituting z = ay/F and v k = zD\M P /\/2, we get 



z" 



vl + [k 2 -—jv x = 0, (48) 

where, prime ' is derivative with respect to conformal time. Summing over all polarization 
states, the Eq. (48) provides us the amplitude of power spectrum of D\ as 

r, t ( 2 W 1 2 2 / H V 1 

P t = 4 x — — 2 -^vl ~ - — -. (49) 



\M 2 J 2ir 2 a 2 F A 7T 2 Vm p ; F' 

So, the ratio of the amplitude of tensor perturbations to scalar perturbations r for /(<3>, i?) 
theories is given by 

r ~ ~ 48e* . (50) 

We now use the measured values of these CMB anisotropy parameters to get the numerical 
values for the parameters (a, b,£, A). 



2.4. Comparison with data 

The recent Planck+WP data [50] gives the curvature perturbation = 2.195to'Hl x 
10" 9 , spectral index n n = 0.9603±0.0073 and r < 0.11(95%CL). We calculate the amplitude 
of the perturbation of the modes which leave the horizon with Nj ~ 70. To get 70 e- 
foldings, the scalar field should roll from 1.6M p to 1M P . Using Eq. (33) and Eq.(36), we 
calculate numerically the values of a and b for different values of self coupling A which 
give the Planck+WP central values of A 2 ^ = 2.195±q;585 x 10 -9 and spectral index n n = 
0.9603 ± 0.0073. We display these correlated set of parameters in Table (I). The slow roll 
parameters are found to be e\ ~ —63 ~ 0.0066 and €2 — —£4 ~ —0.0133 for each value of 
A. The value for the tensor to scalar ratio is also independent of A and turns out to be 
r ~ 0.002. 

The Hubble parameter H in the Jordan frame turns out to be in the range ~ (10~ 5 — 
10~ 4 )M p as A varies in the range (10~ 5 — 0.1). The values of F are found to be much 
larger than unity and hence our assumption of dropping the unity in the expression for F 

is justified. This sets the order of the term ~ 10~ 12 much smaller than e\. 

We have checked that in the limit a ~ 2 and b ~ 1 the correct value of A^ and n-n are 
obtained for A ~ 10" 2 only for large values £ ~ 10 4 . Our Jordan frame calculation in this 
limit is consistent with the results of [25-28] who do the calculation in the Einstein frame. 

We shall next compare the results of these parameters in the Einstein frame and discuss 
the implications of these values in the concluding section. 
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A 


0.1 


io- 2 


10~ 3 


10~ 4 


10~ 5 


a 


3.56398962 


3.27512993 


3.02576840 


2.80955860 


2.62084430 


b 


0.21800513 


0.36243482 


0.48711506 


0.5952182 


0.6895696 


F 


17130 


7286 


2779 


1008 


355 


H 


6.24 x 10~ 4 M P 


4.07 x 10~ 4 M P 


2.5 x 10~ 4 M P 


1.5 x 10~ 4 M P 


8.89 x 10~ 5 Mp 



TABLE I: The values of parameters in the Jordan frame at Nj(<f>j) ~ 70 with £ = 1 and <j)j = 1.6M P 
for different values of A. 



CALCULATION IN THE EINSTEIN FRAME 



Starting with the considered action 



d A X\J—g 



M 2 ( ^R 1 " 1 
" 1 L+ M«+ 2b - 2 



1 A$ 4 

2 M 4 J 



(51) 



we perform a conformal transformation of the metric g^ v to the Einstein frame metric 
g^which is defined as 



~g llv {x) = Q 2 {x)g^(x) , 



where 



n 2 = i + 



£$ a R 



a pfc— 1 



Mp a+2b-2 ■ 



(52) 



(53) 



The Ricci scalar transform as 

R = Q 2 



(54) 



For quasi de-Sitter space we can ignore the second and third terms in the bracket in Eq. 
(54) which is justified in the Eq. (67). For this slow roll case, we can write Eq. (53) in 
Einstein frame as 



n 2 = i + 



m: 



a+2/3 



(55) 



where, a = a/ (2 — b) and f3 = (b — l)/(2 — b). Now we write the action (51) in term of new 
field Xi which is related to the field $ by the relation 



dx 



1 



2£/2 



3o^ 



2q-2" 



1/2 



(56) 
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where £' = £ + p(R/M 2 } Y . This leads the action in term of x as follows 



S f 



d x 



M 2 1 ^ 



where 



For $ >> Mp/£ /1//ct , Eq. (56) can be integrated to give 



(57) 



(58) 



cxp 



2 X 



3 Mp« 



(59) 



Considering g^,, = (— M 2 (t), a 2 (t), a 2 (t), a 2 (t) and varying the action (57) with respect to 
M(t) or a(t) and setting M — 1 in the final equation which corresponds FRW metric, we 
get the Friedmann equation 



12H Z - C _1 M*A I 1 + 



2(3 



a 



0, 



(60) 



where 



( = 12 4/3/c 



W) l /2 (q-2)x 



(61) 



Here we have neglected all the derivative terms of Hubble parameter H . This corresponds to 
slow roll condition, i.e., x 2 is much smaller than potential term. We may write the Hubble 
parameter from Eq. (60) as 



H = M n 



Now using Eq. (62) and (58) we obtain 



[1 + 20/a) A]W« 

— 2(1+13) 
12 £ 



cxp 



2-a \ x 
3 \a + A(3j m 



4/3 



. . 1 ,, 4 , a , 4(l+/3) / 2/3 \ ^ 



cxp 



2l /? /2 - a 



3 V« + 4/3/ M p 



X 



(62) 



(63) 



Here we have taken the approximation exp(^/|-^) >> 1 for x >> M p . We now compute 

the spectral index and curvature perturbation using above potential (63). The slow roll 
parameters for large x » M p comes out to be 



M 2 (U'V 4fa + 2b-4\ 2 - fU" 



I 2 



a + 2b - 4 
3 V a + 46 - 4 



(64) 
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A = 0.1 


A = 10- 2 


A = 10- 3 


A = 10" 4 


A = 10" 5 


a = 3.395 


a = 3.040 


a = 2.750 


a = 2.512 


a = 2.310 


b = 0.2727 


b = 0.4327 


b = 0.5630 


b = 0.6708 


b = 0.7616 



TABLE II: The values of parameters in the Einstein frame at Ne(4>e) — 70 with £ = 1 for different 
values of A. 



and the curvature perturbation 
3H 3 



An = 



2ttU'(x) 



x + 2y + 4 



v 



where x = a + 46 — 4 and y = a + 26 — 4. 
The number of e-folding is calculated as 



1 /3M /xo-Xe 



2 v 2 Vyy V m, 



(66) 



For xo ~ 15M P and Xe ~ 1M P , the number of e-folding is found to be around 70. The 
slow roll parameters e, i] and Hubble parameter H are nearly independent with A and are 
~ 0.02, ~ 0.04 and ~ 5.9 x 10" 5 M p respectively. 

Now from Eqs. (53) and (59), we can calculate the order of terms like Q/Q and (Q/Q) 2 
for » For A = 10" 3 and f = 1, 

fi a(a-l) fA\ / J/^y 

A \4>) 6 \3aM p H J p 

|) , ~^(i)' = 5 .l.xlO-»M; (67) 

respectively, whereas the value of curvature scalar i? = 12if 2 at the same values of parameter 
is 4.24 x 10 _8 M 2 . Here we have neglected <fi as it is small. Thus our approximation made is 
consistent and may be checked for other values of a and b. 



4. RENORMALIZATION 

During inflation the Higgs field is in the Planck regime so the running of the dimensionless 
parameters A and £ due to renormalization from the low energies to the Planck scale must 
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be taken into account [57-64]. What we find in this paper is that in order to reproduce the 
observed density perturbations we need not tune £ to any large values as in ref. [21-28] but 
what is needed is to induce large anomalous dimensions in the scalar and graviton fields [65] 
at the Planck scale. The renoramlization group corrected effective action can be written in 
the general form as 

where t = In with /i being the renormalization scale and the wave-function renormal- 
ization factor 

G * (t) =^-/i^4 (69) 

where 7$(t) is the anomalous dimension of the scalar field. 

A similar expression holds for the graviton wave- function renormalization G g . The pa- 
rameters a and b can formally be related to wave-function renormalizations G$ and G g 

as 

n dlnGq, , „ n dlnG a 

a = 2 + 2 w ; " = 1 + 2 mC- (70) 

In this paper we show that the correct values of a and b which generate a viable inflation 
must lie in the range a = (2.3 — 3.6), b = (0.68 — 0.22) and we believe that renormalization 
effects should lead to these values of anomalous dimension of the scalar and graviton fields. 



5. DISCUSSION AND CONCLUSION 

We have generalised the curvature coupling models of Higgs inflation to study inflation 
with a scalar field for a |<3> 4 potential and a curvature coupling of the form . This 

curvature coupling can be obtained for example by starting with a tree level coupling term 
£$ 2 i? (a ~ 2, b ~ 1) if the wave-function renormalization of $ and the graviton introduce 
large anomalous dimensions for these fields such that a is increased to lie in the range 
(2.3 — 3.6) and b is reduced to the range (0.68 — 0.22). In this case if the renormalized A lies 
in the range (10~ 5 — 0.1) at the Planck scale, one can have a viable inflation model even for 

It has been shown that the Higgs self coupling can go from A = 0.13 at the electroweak 
scale for the 125 GeV Higgs to A = 10 -5 at the Planck scale by tuning the top mass or by 
introducing extra interactions [15, 16, 49]. This leads us to conclude that the Higgs field 
may be still a good candidate for being the inflaton in the early universe if one considers 
higher order curvature coupling which may arise from renormalization at the Planck scale. 

The tensor to scalar ratio r in this model is small and the with scalar curvature 
couplings is not ruled out by observational limits on r unlike the pure |$ 4 theory [8, 66]. 
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We find that the results of Jordan frame and Einstein frame calculations are in agreement 
up to leading order in the slow roll parameters as expected from the approximation made 
in the course of conformal transformation (54,55) (and justified in the Eq. (67)). Finally, 
by requiring the curvature coupling parameter to be of order unity, we have evaded the 
problem of unitarity violation in scalar-graviton scatterings [29] which plagued the £$ 2 i? 
Higgs inflation models [25-28]. 
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